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The properties of the new analytic running coupling are investigated at the higher loop 
levels. The expression for this invariant charge, independent of the normalization point, 
is obtained by invoking the asymptotic freedom condition. It is shown that at any loop 
level the relevant {3 function has the universal behaviors at small and large values of the 
invariant charge. Due to this feature the new analytic running coupling possesses the 
universal asymptotics both in the ultraviolet and infrared regions irrespective of the loop 
level. The consistency of the model considered with the general definition of the QCD 
invariant charge is shown. 

1. Introduction 

The description of hadron interaction in the infrared region remains an actual 
problem of Quantum Chromodynamics (QCD). Since the standard perturbation 
theory provides no definite answer on this problem, a variety of nonperturbative 
methods is usually invoked for the comprehensive investigation of this matter. The 
point is that the renormalization group (RG) summation leads to the violation of the 
proper analytic properties of the relevant physical quantities, that contradicts the 
general principles of the theory. For instance, the ghost pole appears in the expres- 
sion for the running coupling at the one-loop level. Therefore, in order to improve 
this situation one has to involve into consideration the condition of analyticity. 

Being based on the first principles of the local Quantum Field Theory, the an- 
alytic approach seeks to recover the violated after RG summation proper analytic 
properties of the relevant physical quantities. Its original ideas were formulatedtra 
in the framework of Quantum Electrodynamics (QED) in the late 1950's. The ba- 
sic idea of this approach is an explicit imposition of the causality condition which 
implies the requirement of the analyticity in the q 2 variable for the relevant .phys- 
ical quantities. The analytic approach has recently been extended to QCDu and 
applied to the 'analytization' of the perturbative series for the QCD observables.Q 
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The term 'analytization' means the recovering of the proper analytic properties in 
the q 2 variable by making use of the Kallcn Lchmann spectral representation 
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with the spectral density g(a) determined by the initial (perturbative) expression 
for a quantity A(q 2 ): 

q(o) = — lim [a(-<t - is) - A(-a + is)] , a>0. (2) 
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The analytic approach has recently been applied to the 'improvement' of the 
(3 function perturbative expansion .oB In accordance with the model proposed in 
Refs. U,0, the new analytic running coupling (NARC) is the solution of the renor- 
malization group equation for the invariant charge with the anajvtized (3 function. 
This running coupling possesses a number of appealing features.Q'Q Its most impor- 
tant advantages are the following. First of all, the new analytic running coupling has 
no unphysical singularities. Further, the NARC explicitly involves the ultraviolet 
(UV) asymptotic freedom with the infrared (IR) enhancement in a single expres- 
sion. It is worth noting here that there is a number of evidences for such a behavior 
of the QCD invariant charge. In particular, the recent lattice simulationsfE as well 
as the solution of the Schwinger-Dyson equationsEJ testify to the IR enhancement 
of the QCD running coupling. Furthermore, it was demonstrated recently&El that 
the NARC provides the quark confinement, the one-gluon exchange model being 
employed. Remarkably, the new model for the analytic running coupling has no 
additional parameters, i.e., similarly to the perturbative approach, Aqcd remains 
the basic characterizing parameter of the theory. 

In this letter the investigations of the new analytic running coupling is contin- 
ued. The primary objective is to study the properties of the NARC at the higher 
loop levels. Since we have only the integral representation for the NARC here, its 
straightforward investigation turns out to be rather complicated. Nevertheless, by 
examining the properties of the (3 function corresponding to the NARC one suc- 
ceeded in the description of the asymptotic behavior of the new analytic running 
coupling at the higher loop levels. 

The layout of the letter is as follows. In Sec. 2 the new analytic running coupling 
is briefly discussed. By invoking the asymptotic freedom condition the expression 
for the NARC, independent of the normalization point, is obtained. Further, the 
conclusions concerning the higher loop and scheme stability of the current approach 
are drawn. In Sec. 3 the [3 function corresponding to the new analytic running 
coupling at the higher loop levels is investigated. It is shown that this j3 function 
has universal behaviors both at the small and large values of the invariant charge, 
irrespective of the loop level considered. This results in the universal asymptotics 
of the new analytic running coupling both in the UV and IR regions at any loop 
level. In Sec. 4 the compatibility of the new model for the analytic running coupling 
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with the general definition of the QCD invariant charge is shown. The ways of the 
analyticity requirement implementation are discussed. In Sec. 5 the obtained results 
are summarized. 



2. New Analytic Running Coupling in QCD 

The complementation of the (3 function perturbative expansion with the analyt- 
icity requirement leads to the analytized RG equation for the new analytic running 
coupling (see Refs. |[0 for the details). At the £-loop level this equation acquires 
the form 

dhv — --sz^H^ ^n ' B > -^+T- (3) 

\-* /an 

Here N ain(/i 2 ) is the new analytic running coupling at the Moop level, ai\^i 2 ) 
denotes the £-loop perturbative running coupling, 3(/i 2 ) = a(/i 2 ) /3 /(47r), f3 = 
11 — 2 rif 1 3 and (3% = 102 — 38 n{ / 3 are the standard coefficients for the (3 function 
expansion, and nf is the number of active quarks. At the one-loop level Eq. (|§|) can 
be integrated explicitly with the result! 

^V) = ?^. * = (4) 



The properties of the one-loop NARC (Q) have been studied in details in Ref. ||. 
It is worth to emphasize here that the new analytic running coupling (|i|) possesses 
the following significant feature 

N «i 1 n ) (9 2 ) = f N «^(^)- (5) 

Recently it has been shownEI that this symmetry of the QCD invariant charge 
precisely corresponds to the conformal inversion symmetry of the instanton size 
distribution.^! In particular, this implies that one can also derive the expression (|J) 
for the QCD running coupling proceeding from the entirely different motivations. It 
is interesting to note also that the relation 
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holds for all q 2 > in the framework of the model proposed. Here N ain (q 2 ) is the 
invariant charge (||) and N /?in denotes the relevant (3 function (see Eq. (|l^) further). 

Let us turn to the higher loop levels. From the very beginning one should 
note that the solution of Eq. (||) is determined up to a constant factor due to the 
logarithmic derivative on its left-hand side. In previous studies 3, this problem has 

a In Ref. the factor zq/z has been taken out from the exponent of Eq. (Jt|) due to the property 

P° TL^Ha) a" 1 da = 1 for all I. 
J 
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been eliminated by normalization of the solution of Eq. (^) to its value at a point c; 2 

da 
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However, in some cases it turns out to be more convenient to deal with the explicit 
expression for the running coupling independent of the normalization point. 

In this letter we propose a simple physical method for removing the ambiguity 
mentioned above. Indeed, this can easily be achieved by involving the condition 
of the asymptotic freedom, b namely ™a&2(q 2 ) — > ai\q 2 ) when q 2 — > oo (in fact, 
this has already been used in Eq. ([|)). It is worth noting that this method does 
not violate the renormalization invariance of Eq. (0). In general, one becomes 
able to derive the integral representation for the Moop new analytic running cou- 
pling, independent of the normalization point, by invoking the similar condition 
N 5iri(<Z 2 ) — > (<Z 2 ) when q 2 — > oo. For this purpose let us consider the ratio 
of the £-loop new analytic running coupling (Q) normalized at a point qfi to the 
one-loop NARC written in the form (fy and normalized at the same point q 2 : 
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where ATZ w (a) = TZ {l) (a) - TZ {1) {a). Proceeding to the limit q 2 
at the expression for the £-loop new analytic running coupling 



oo, we arrive 
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It is worth noting that there is the integral representation of the Kallen-Lehmann 
type for the new analytic running coupling 
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where 
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In particular, this was used in Ref. ^ when evaluating the normalization coefficients. 
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is the £-loop spectral density, 



V cr J In a + -k z 

is the one-loop spectral density, and 

1^V)=7r pXRViOf. (14) 

In the exponent of Eq. (|l2|) the principle value of the integral is assumed. However, 
in the practical use the expression ([lO]) turns out to be more convenient (at least, 
beyond the one-loop level). 

One of the important features of the new analytic running coupling is the absence 
of unphysical singularities at any loop level. Furthermore, it involves both the 
asymptotic freedom behavior and the IR enhancement in a single expression. This 
feature enables one to describe a broad range of physical phenomena including both 
the standard perturbative and the intrinsicallynonperturbative ones (see Ref. || for 
the details). In particular, it has been shownErl!3 that in the framework of the one- 
gluon exchange model the new analytic running coupling (Q) explicitly leads to the 
rising at large distances static quark- antiquark potential. Let us emphasize here that 
the additional parameters are not introduced in the theory. The detailed description 
of the properties and advantages of the NARC (0) is given in the papers.Qa 




Fig. 1. The new analytic running coupling at different loop levels, z = q 2 /A 2 . 

The Figure |l| shows the new analytic running coupling ( |Io|) at the one-, two- 
and three-loop levels. It is obvious that NARC possesses the higher loop stability. 
Thus, the curves corresponding to the two-, and three-loop levels are practically 
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indistinguishable. Proceeding from this one can also draw the conclusion concerning 
the scheme stability of the current approach. In particular, in Fig. pi the curve 
corresponding to the three-loop new analytic running coupling N a4 n is platted by 
making use of the coefficient 02 — 2857/2 — 5033 rif/18+325 /54 computedtj in the 
MS scheme. The account of the third term on the right-hand side of Eq. (|^) does not 
lead to a valuable quantitative variation of its solution in comparison with the two- 
loop approximation. Therefore it is clear that using the coefficient computed in 
another subtraction scheme c does not lead to significant variation of the solution to 
Eq. (||) in comparison with the considered case of the MS scheme. This statement 
follows also from the fact that the contribution of every subsequent term on the 
right-hand side of Eq. (||) is substantially suppressed by the contributions of the 
preceding ones. 

Since we have only the integral representation for the new analytic running 
coupling at the higher loop levels, its straightforward investigation becomes rather 
complicated. Nevertheless, the examining of the (3 function corresponding to the 
NARC enables one to elucidate a number of important questions, in particular, the 
asymptotic behavior of the new analytic running coupling. 



3. The (3 Function: Higher Loop Levels 

In the previous letter! the function 

d lna(> 2 ) 

p[a) = (15) 

corresponding to the one-loop new analytic running coupling (|]) has been derived d 
and studied in details. Due to the explicit expression for the NARC at the one- loop 
level, one succeeded in performing the relevant investigation manifestly. Thus, the 
(3 function (|l5|) corresponding to the one-loop NARC (^) acquires the form 

1 - N(a)/a 



where the function N(a) 
N(a)- 



ln[N(a)} ' 
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iV-i(a), 1 < a, 
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is defined in terms of the many-valued Lambert W function' 

W k (x) exp [Wk (a;)] = x 



exp 



(16) 



(17) 



(18) 



and the result 
(q 2 ) = 1/lnz, 



(see Ref. |8J for the details). Figure |2j presents the (3 function (|l 
corresponding to the one-loop perturbative running coupling a 
namely (3^ {a) = —a. 

c At least, in schemes that do not have unnaturally large expansion coefficients (see Ref. ^ and 
references therein for the detailed discussion of this matter), 
d In Ref. ^the definition flip ) = da(fi 2 )/d ln^t 2 was used. 

e In definitions (|17|) and (Jl 8[) k denotes the branch index of the Lambert W function. 
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Fig. 2. The f3 function corresponding to the one-loop new analytic running coupling (solid curve). 
The relevant perturbative result is shown by the dot-dashed line. 



The investigation!] of the properties of the function N(a) defined in Eq. (|Ti 
enables one to find the explicit asymptotic behavior of the j3 function (|T6|). Thus, 
for small values of the running coupling ^ain Eq. (|l^) coincides with the well-known 
perturbative result 



-a + O 



a exp 
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(19) 



It is worth noting here that the second term in Eq. (|19|) points to the intrinsically 
nonpcrturbative nature of the (3 function (hq) . For large values of a one has 



-1 + 



ln(lna) 
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Such a behavior of the (3 function provides the IR enhancement of the invariant 
charge, namely (up to a logarithmic factor) a(q 2 ) ~ 1/z, when z — ► 0. In turn, this 
leads in a straightforward way to the confining static quark-antiquark potential 
(see Refs. |^, [l2]for the details). Therefore, at the one- loop level the (3 function ( |l6| ) 
explicitly incorporates the UV asymptotic freedom and the IR enhancement of the 
running coupling in a single expression. 

Let us proceed to the higher loop levels. As it has been mentioned in the previous 
section, here we have only the integral representation for the new analytic running 
coupling ([n]) . This fact significantly complicates the investigation and leads to the 
necessity of applying the numerical methods. It is worth to mention here again 
that we originate in the standard perturbative expansion for the [3 function and 
complement it with the analyticity requirement. The results of the computations 
corresponding to the two-, and three-loop levels are shown in Figs. [| and |], respec- 
tively. Figure || presents the (3 function N /?i n ' ) (a) corresponding to the two-loop new 
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I ■ : ■ . . 3. The /3 function corresponding to the two-loop new analytic running coupling (solid curve). 
The relevant perturbative result is shown as the dot-dashed curve. 

analytic running coupling N a^(q 2 ) together with the respective perturbative result 
/3s (a) = —a — B\a 2 . Figure [| shows the analogous functions at the three-loop 
level, namely N /?iJ (a) and 0s (a) = —a — Bid 2 — B 2 a 3 . It is clear from the figures 
||, [|, and ^ that the f3 function corresponding to the new analytic running coupling 
coincides with its perturbative analog in the region of small values of the invariant 
charge at any loop level. In other words, in the framework of the model in hand 
the complete recovering of the perturbative limit in the UV region takes place. 




Fig. 4. The function corresponding to the three-loop new analytic running coupling (solid curve). 
The relevant perturbative result is shown as the dot-dashed curve. 

Let us turn now to the asymptotics of the (3 function N /?in (a) at large values 
of the running coupling. Figure |s| presents the /3 functions N /3in(a), £ = 1,2,3 



8 



0.0 



-0.5 



-1.0 



-1 .5 



Vpj(a 

• 


) 

2,3 loops 


| l 


l0 °P log 10 a 



10 12 



Fig. 5. The j3 function corresponding to the new analytic running coupling at different loop levels 
(solid curves). The one-loop perturbative result is shown as the dot-dashed curve. 



and the one-loop perturbative result /3s (a) = —a. This figure clearly shows the 
perturbative limit at small a, as well as the universal asymptotic behavior of the (3 
function corresponding to the NARC N pin(a) — ► — 1, when a — > oo. 

The latter statement can also be proved in an independent way. Indeed, due 
to the IR enhancement of the new analytic running coupling, the value of the 
right-hand side of Eq. (j|) (it is nothing but the f3 function at the relevant loop 

level) when fi 2 — > corresponds to the limit N ain — ► oo. One can show that 

i+i 



lim , 



^ J a 



= 1, and lim. 



{ ai e \q 2 ) | = (j > 1 is a positive 
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integer number) , irrespective of the loop level. Therefore, we infer that at any loop 
level the (3 function corresponding to the new analytic running coupling tends to 
the universal limit 

lim %^{a) = -l. (21) 

a — >oo 

As it has been mentioned above, such a behavior of the (3 function leads to the IR 
enhancement of the invariant charge, namely a(q 2 ) ~ 1/z, when z ^ 0. 

Thus, the f3 function corresponding to the new analytic running coupling has 
the universal asymptotic behaviors both at the small ( N /3in(a) ~ —a) and large 
( N f3in(a) ~ — 1) values of the invariant charge irrespective of the loop level. There- 
fore, the new analytic running coupling possesses the universal asymptotics both 
in the UV ( N a ( ^(q 2 ) ~ 1/hxz) and IR ( N a { ^{q 2 ) - 1/z) regions at any loop level. 
In particular, this implies that using the new analytic running coupling (|Io| ) at the 
higher loop levels will also results in the confining quark-antiquark potential. 
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4. Discussion 

In general, the QCD invariant charge is defined as the product of the corre- 
sponding Green functions and vertexes. For example, in the transverse gauge the 
following definition 

a(q 2 )=a( t x 2 )G(q 2 )g 2 (q 2 ) (22) 

takes place, where a{q 2 ) is the QCD running coupling, a{ji 2 ) denotes its value 
at a normalization point /i 2 , G{q 2 ) and g{q 2 ) are dimensionless gluon and ghost 
propagators, respectively. It is essential to note here that the Green functions 
possess the proper analytic properties in the q 2 variable (namely, there is the only 
left cut q 2 < 0) before the RG summation. Indeed, for the considering in this 
letter massless case a Green function can be represented as the perturbative power 
series in a{y?) with the coefficients being a polynomial in ln(g 2 /^ 2 ). Obviously, 
in this case the integral representation of the Kallen-Lehmann type must hold for 
any finite order of the perturbative expansion for the Green function/ The RG 
summation leads to violation of such a properties of the Green functions (e.g., 
in the simplest one-loop case the ghost pole appears). As it has already been 
mentioned above, the analytic approach to QCD seeks to recover the proper analytic 
properties of the relevant quantities. However, for the consistency of involving the 
analyticity condition with the definition of the invariant charge (p2|), one has to 
apply the analytization procedure to the logarithm of the Green function. Really, if 
the spectral function of the Kallcn Lchmann representation for the Green function 
G is of a fixed sign (that is true for the leading orders of perturbation theory), then 
G has no nulls in the complex q 2 plane. Hence, its logarithm, In G, can also be 
represented as the spectral integral of the Kallen-Lehmann type. 

In the framework of perturbative approach the considered above Green func- 
tions have the following form: g(q 2 ) = (1/ln z) d <> and G(q 2 ) = (1/ln z) d ° , where 
z = q 2 /A 2 , d g and cLq denote the corresponding anomalous dimensions. For the 
latter the relationship 2d g + do = 1 holds, that plays the crucial role in the defini- 
tion ( p2| ) . Applying the analytization procedure to these functions in the described 
above way, one arrives at the following result: g an (q 2 ) = [(z — l)/(z ln2;)] ds and 
G an {q 2 ) = [(z - l)/(z \nz)] dG . Therefore, in this case the invoking the analyticity 
condition does not affect the definition of the invariant charge (p2|). It is inter- 
esting to note here tha.t the similar situation takes place in the nonperturbative 
a-expansion methodllj'EZl also. Thus, we infer that the NARC ^ is consistent with 
the general definition of the QCD invariant charge. 

It is worth to mention that there are different methods of removing of the un- 
physical singularities from the running coupling in the framework of the analytic 
approach to QCD. In original workcl the analytization procedure was applied to 
the perturbative invariant charge in a straightforward way (see Ref. |falso). The 
model in handES is based on the complementation of the (3 function perturbative 

' One should note that, unlike the case of the QED photon propagator, the relevant spectral 
function here can't be directly identified with the spectral density. 
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expansion with the analyticity condition. In both cases the additional parameters 
are not introduced into the theory, the models being the 'minimal' ones in this sense. 
Of course, in the UV region these models have identical behavior determined by 
the asymptotic freedom. However, there is a qualitative distinction between them 
in the IR region (see discussion in RefjT] also). 

Let us note here that in equation (0) for the new analytic running coupling the 
analytization of its right-hand side as a whole is assumed. In turn, this results in 
the representation of t he co rresponding (3 function as a non-power series. As it 
has already been shown,@@ it is this way of the analyticity requirement involving 
that leads to the stability of obtaining results in a broad range of energies. In 
particular, such an implementation of the analytization procedure ultimately leads 
to the universal IR behavior of the NARC. 

As it was mentioned above, the NARC incorporates the IR enhancement with 
the UV asymptotic freedom in a single expresskm. The recent results of the nonpiir- 
turbative studies (namely, lattice simulationsQ't3 and Schwinger-Dyson equationstll) 
testify to such a behavior of the QCD invariant charge. For the completeness of the 
pattern let us also note here that these nonperturbative methods, being the matter 
of contemporary comprehensive investigations, provide no unique point of view on 
the IR behavior of the QCD invariant charge. 

5. Conclusions 

The properties of the new analytic running coupling are studied at the higher 
loop levels. By making use of the asymptotic freedom condition the expression for 
the NARC, independent of the normalization point, is obtained. The conclusions 
about the loop and scheme stability of the current approach are drawn. The (3 
function corresponding to the NARC is constructed and examined. It is shown that 
the behaviors of this (3 function at both its key asymptotics (namely, when a — » 
and when a — ► oo) are the same at any loop level. This results, irrespective of 
the loop level, in an universal asymptotics of the new analytic running coupling. 
Namely, its UV behavior ( N 5in(g 2 ) ~ 1/lnz, when q 2 — > oo) is determined by the 
asymptotic freedom, and there is the IR enhancement of the new analytic running 
coupling ( N 5&2 (q 2 ) ~ 1/z, when q 2 — > 0). The consistency of the model considered 
with the general definition of the QCD invariant charge is shown. 

Acknowledgments 

A. N. would like to thank Prof. F. Schrempp for the stimulating comments. 
I. S. is grateful to Profs. D. I. Kazakov and D. V. Shirkov for interest in this work. 
The partial support of RFBR (grants 99-01-00091, 99-02-17727, and 00-15-96691) 
is appreciated. 



11 



References 



11. 
12. 

13. 
14. 
15. 
16. 
17. 
IS. 
19. 
20. 



hep-ph/9704333 . 
(1999); Sep- 



1220 



P. J. Redmond, Phys. Rev. 112, 1404 (1958). 

N. N. Bogoliubov, A. A. Logunov, and D. V. Shirkov, Zh. Eksp. Tear. Fiz. 37, 805 
(1959) [Sov. Phys. JETP 37, 574 (I960)]. 

D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79, 1209 (1997) 
I. L. Solov tsov and D. V. Shirkov, Theor. Math. Phys. 120, 
ph/990930s| . 

Phys. Rev. D62, 094028 (2000); |hep-ph/9912351 . 

in Proceedings of the Sixth Workshop o n Non-Perturbati ve Quantum 
Paris, France, 2001 (to be publi shed) ; 



A. V. Nesterenko 
A. V. Nesterenko. 
Chromodynamics, 
A. V. Nesterenko 
A. V. Nesterenko 



hep-ph/0106305 



Phys. Rev. D64, 116009 (2001); |hep-ph/0102124|. 
Mod. Phys. Lett. A15, 2401 (2000); |hcp-ph/0102203| . 
UKQCD Collabo ration, D. A. Smith and M. J. Teper, Phys. Rev. D58, 014505 (1998); 
hep-lat/98010081 . 



10. ALPHA Collabo ration, A. Bode et ai, Phys. Lett. B515, 49 (2001), hep-lat/0105003 



J. Heitger et ai, |hep-lat/011020l| . 

A. I. Alekseev and B. A. Arbuzov, Mod. Phys. Lett. A13, 1747 (1998); |hep"-ph/970422£ . 
A. V. Nesterenko, in Proceedings of the 4th International Conferenc e on Quark Con- 
finement and the Hadron Spect rum, Vienna, Austria, 2000 (in press); |iep-ph/0010257 . 
F. Schrempp, |hep-ph/0109032| . 

0. V. Tarasov, A. A. Vladimirov and A. Yu. Zharkov, Ph ys. Lett. B93, 429 (1980). 
P. A. Raczka, Z. Phys. C65, 481 (1995); |hep-ph/940734; . 

1. L. Solovtsov, Phys. Lett. B327, 335 (1994); B340, 245 (1994). 

A. N. Sisakyan and I. L. Solovtsov, Phys. Part and Nucl. 30, 4 61 (1999). 

D. V. Shirkov, Eur. Phys. J. C22, 331 (2001); |hep-ph/0 10728^ . 

I. L. Solovtsov and D. V. Shirkov, Phys. Lett. B442, 344 (1998); |hcp-ph/9711251 . 

K. A. Milto n, I. L. Solovtsov , O. P. Solovtsova, and V. I. Yasnov, Eur. Phys. J. C14, 

495 (2000); ^iep-ph/0003030] . 



12 



